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In the time evolution of neutral meson systems, a perfect quantum coherence is usually assumed. 

The important quantities of the B° system, such as sin 2/3 and Amj, are determined under this 
assumption. However, the meson system interacts with its environment. This interaction can 
lead to decoherence in the mesons even before they decay. In our formalism this decoherence is 
modelled by a single parameter A. It is desirable to re-examine the procedures of determination of 
sin 2/3 and A uid in meson systems with decoherence. We find that the present values of these two 
quantities are modulated by A. Re-analysis of B° data from B-factories and LHCb can lead to a 
clean determination of A, sin 2 (3 and A mu- 


I. INTRODUCTION 

In neutral meson systems, quantum coherence plays 
a crucial role in the determination of many observables. 
However, any real system interacts with its environment 
and this interaction can lead to a loss of quantum coher¬ 
ence. The environmental effects may arise at a funda¬ 
mental level, such as the fluctuations in a quantum grav¬ 
ity space-time background Hi- They may also arise 
due to the detector environment itself. Irrespective of 
the origin of the environment, its effect on the neutral 
meson systems can be taken into account by using the 
ideas of open quantum systems M- This formalism 
enables the inclusion of effects such as decoherence and 
dissipation in a systematic manner . Such an inclusion 
is in accordance with the general principle of fluctuation- 
dissipation theorem which states that dissipation is bal¬ 
anced by fluctuations. 

The time evolution of neutral mesons, which are coher¬ 
ently produced in meson factories, are used to measure 
a number of parameters of the standard model of par¬ 
ticle physics and also to search for physics beyond the 
standard model. However, decoherence is an unavoid¬ 
able phenomenon as any physical system is inherently 
open due to its inescapable interactions with a pervasive 
environment. With the inclusion of the decoherence ef¬ 
fects, the measured values of some of these parameters 
can get masked. As the source of decoherence in the case 
of mesons could be expected to be coming from a much 
finer scale, it may happen that the numerical value of 
some of the masked observables are not greatly affected. 
This should however be verified experimentally. 

In this work, we study the effect of decoherence on 
the important observables in the B° meson system, such 
as the CP violating parameter sin 2/3 and the 5° — B ° 
mixing parameter Amj. We show that these parameters 
are affected by decoherence. So far only one attempt 
has been made to determine decoherence in Bd meson 
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system Q. The bounds on the decoherence parameter 
were obtained from the data on Rd, the ratio of the total 
same-sign to opposite sign dilepton rates in the decays of 
coherent Bd — Bd coming from the T(4S * I ) decays. The 
data on Rd has not been updated in the last two decades 
[U, whereas the B-factories have provided direct and pre¬ 
cise information on the Bd — Bd mixing parameters. In 
this work, we also suggest a number of methods which 
will enable clean determination of the decoherence pa¬ 
rameter along with the other observables quite easily at 
the LHCb or B-factories. We also attempt determination 
of the decoherence parameter and A rrid using Belle data 
on the time dependent flavor asymmetry of semi-leptonic 
Bj decays as given in Ref. [§j. 

The evolution of the B° system is built up from first 
principles. The effect of the environment forces the evo¬ 
lution to be a semi-group rather than a unitary one 
10, EH- We use the density matrix formalism to rep¬ 
resent the time evolution of the B® system. This ensures 
the complete positivity of the state of the system and 
hence its physical validity. In this formalism, the deco¬ 
herence is modelled by a single parameter A. By con¬ 
struction, the density matrices are trace preserving. 

The work presented here, we hope, would lead to the 
inclusion of the effects of decoherence in the analysis of 
data from the B® systems. It may be worthwhile to rean¬ 
alyze the data from the B factories and LHCb to verify 
if a signature of decoherence is already inherent in it. 
Thus a detailed study of B° observables can lead to tests 
of physics at scales much higher than those typical of 
flavour physics. 

We first study the parameter sin 2/3, whose measure¬ 
ment is the first signal for CP violation outside the neu¬ 
tral kaon system. The precision measurement of its value 
is the corner stone in establishing the CKM mechanism 
for CP violation. With the inclusion of the decoherence 
effects, it turns out that the experimentally measured CP 
asymmetry depends both on the decoherence parameter 
A and the angle /? of the unitarity triangle. Next we 
study Arrid, which denotes the mixing in the B® system 
and is an important input in extracting sin 2/3 from the 
measured time dependent CP asymmetry. We find that 
Arrid is also affected by the decoherence effects. Finally, 
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we suggest a method of analysis by which the three quan¬ 
tities, (a) A, (b) A rrid and (c) sin 2/3 can all be measured. 


II. DETERMINATION OF sin 2d 

In the following, we develop the formalism which is ap¬ 
plicable to B ° as well as 13° mesons. We are interested 
in the decays of B° and B° mesons as well as B° o B° 
oscillations. To describe the time evolution of all these 
transitions, we need a basis of three states: |-B°), |l3 0 } 
and |0), where |0) reprents a state with no B meson and 
is required for describing the decays. In this basis, we 
can define Pbo(bo)( 0), the initial density matrix for the 
state which starts out as B°(B°). The time evolution of 
these matrices is governed by the Kraus operators i<Q(f) 
as p(t) = Ai(t)p(O)A'J(f) 0. The Kraus operators 
are constructed taking into account the decoherencc in 
the system which occurs due to the evolution under the 
influence of the environment 00- The time depen¬ 
dent density matrices are 

0 

0 

2(e rt - a ch ) 

0 
0 

2(e rt - a ch ) 

for B u and 13°, respectively. In the above equation, 
a c h = cosh (ATI), ash = sinh(A£l) ; 0c = cos (A mt), 
a s = sin (Amt), T = (T^ + V H )/ 2, AT = T L - 
where T l and T h are the respective decay widths of the 
decay eigenstates B° L and B^. Also A is the decoherence 
parameter, due to the interaction between one-particle 
system and its environment. As our main motivation 
is to bring out the fact that fundamental parameters of 
B — B mixing and B sector CP violation are affected by 
decoherence, here, and from here on, we will neglect the 
small mixing induced CP violation to keep our formulae 
simple. 

We define the decay amplitudes Af = A(B° —> /) and 
Af = A(B° —> /). The hermitian operator describing 
the decays of the B° and B° mesons into / is 

( \A f \ 2 A f *A f 0 \ 

Of = A f A* \Af\ 2 0 . (2) 

\ 0 0 0 / 

The probability, Pf(B°/B°;t), of an initial B°/B° 

decaying into the state / at time t is given by 

Tr [Of PB°(B°)(i)]- 

Let us now consider B d —»■ J/i^Ks decay. One can 
define a CP violating observable 


Pbo ( t) 
ip-rt 

ip-rt 

2 e 


® ch & CL c 

-a 8 h - ie~ xt > 

&sh a s 

a c h - e~ xt a, 

0 

0 

a c h ~ e~ xt a c 

—o>sh + ie~ Xt ' 

® sh ^ 0'S 

a c h + e~ xt a, 

0 

0 


Aj/^Ks (t) 


Pj/^Ks (P° ;t)- Pj/ipK s {B Q d ; t ) 
Pj/^K s (Pd'A) + Pj/^Ks(B° d - : t) 


(3) 


Calculating the probabilities using Eqs. (ID and © we 
get 

Aj/^Ksi*) _ (| A/1 2 - l) cos(Am d f) + 2Im(A/) sin (Am d f) 
^ ~ (1 + |A/1 2 ) cosh (AAA) _ 2Re(A/) sinh (AAil) ’ 

(4) 

where A/ = A{B° d ->• J/%l)K s ))IA(B 0 d ->■ J/ipK s ). 
Putting A = 0 in the above equation, we get the usual ex¬ 
pression for CP asymmetry in the interference of mixing 
and decay. Thus the presence of decoherence modifies 
the expression for CP asymmetry in the interference of 
mixing and decay. 

In order to determine sin 2/3 from asymmetry defined 
in Eq. 01 it is usually assumed that, Ar d 0, |A/| = 1, 
i.e., no direct CP asymmetry and Im(A/) ~ sin 2/3. With 
these approximations, the above expression simplifies to 

Aj/^Ks W = sin 2 /3 e -At sin (A m d t) . (5) 

Therefore we see that the coefficient of sin (Am^t) in 
the CP asymmetry is sin2/3e^ At and not sin 2/3! The 
measurement of sin 2/3 is masked by the presence of de¬ 
coherence. Thus in order to have a clean determination 
of sin 2/3, an understanding of A is imperative. 

I , Decoherence is expected to come from a scale much 
finer than that of flavor physics and is likely to be small. 
Therefore, in the actual comparison to the data, one 
Lshould include all the known effects, which are usually 
I neglected in the extraction of sin 2/3 and then do a fit 
for clean determination of sin 2/3 and A. The full fledged 
formula, of course, will include the CP violation in mix¬ 
ing and decay width AT^. Apart from these effects, one 
should also take into account the penguin contributions. 

The theoretical precision for the extraction of CP violat¬ 
ing phase sin 2/3 from the CP asymmetry of i3° —> J/i/jKg 
decay, defined in Eq. 01 is limited by contributions from 
doubly Cabibbo-suppressed penguin topologies 0 , 0 . 
This involves computation of non-perturbative hadronic 
parameters which, at present, cannot be achieved reli¬ 
ably using QCD. However, a way to control the penguin 
effects is offered by the U-spin symmetry of strong in¬ 
teractions which relates B° s —»■ J/ipKs to B d —> J/ipKs 
0. Ref. [16] discusses the constraining of the relevant 
penguin parameters by making use of this symmetry as 
well as plausible assumptions for various modes of similar 
decay dynamics. 


III. DETERMINATION OF A m d 

It is obvious that in order to determine sin 2^, we need 
to know A nid and A. If A nid is measured using observ¬ 
ables which are independent of A, then we only need to 
determine A for the clean extraction of sin 2/3. If the de¬ 
termination of A md is also masked by the presence of 
decoherence then we need to have a clean determination 
of A md- 

The present world average of Aquoted in PDG is 
(0.510 ± 0.003) ps -1 0 which is an average of measure- 
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ments of A m d from OPAL 0, ALEPH 0L DELPHI 
0, L3 0LCDF 0, BaBar 0, Belle lH, DO 0 
and LHCb [271 ] experiments. There are several ways in 
which A md can be determined experimentally. LHCb, 
CDF and DO experiments determine Am^ by measuring 
rates that a state that is pure B° at time t = 0, decays 
as either as B® or B ° as function of proper decay time. 
In the presence of decoherence, the survival (oscillation) 
probability of initial B® meson to decay as at a 

proper decay time t is given by 


e -rt 

P±(t, A) = —— [cosh(APrft/2) ± e At cos(Amdt)] ■ 

( 6 ) 

The positive sign applies when the B® meson decays with 
the same flavor as its production and the negative sign 
when the particle decays with opposite flavor to its pro¬ 
duction. We see that the survival (oscillation) probabil¬ 
ity of B d is A dependent! The time dependent mixing 
asymmetry, used to determine Am^, is then given by 


P+{t, A) - P_(t, A) _ _ xt cos(A m d t) 

mixl ’ ’ P+(i, A) + P~(t, A) ' cosh(Ar d t/2)' 

(7) 

Thus we see that the in the limit of neglecting Ar^, the 
otherwise pure cosine dependence of mixing asymmetry 
is modulated by e~ xt . Belle and BaBar experiments de¬ 
termine Awd by measuring time dependent probability 
P+(t) of observing unoscillated B°B° events and P~{t) of 
observing oscillated B d B d /B^B d events for two neutral 
Bd mesons produced in an entangled state in the decay 
of the T(4S) resonance. The expressions for P±(t), in 
the presence of decoherence, are the same as those given 
in Eq. m, except that the proper time t is replaced by 
the proper decay-time difference At between the decays 
of the two neutral Bd mesons. Therefore, we see that 
the determination of A md at LHCb, CDF, DO, Belle and 
BaBar experiments is masked by the presence of A. The 
true value of A md, along with AT^, can be determined 
by a three parameter (A md, AT^, A) fit to the time de¬ 
pendent mixing asymmetry A m i x (f, A) defined in Eq. (|T|). 
This in turn will enable a determination of true value of 
sin 2/3 using Eq. 6D- 

Determination of A md in the LEP experiments is 
mainly based on time independent measurements, i.e., 
from the ratio of the total same-sign to opposite-sign 
semileptonic rates (Rd) or the total B^ — B® mixing prob¬ 
ability (xd)■ We shall now see that these observables are 
also A dependent. Therefore all the methods used to de¬ 
termine A md depend upon A. 


IV. CORRELATED B° d MESON 
SEMI-LEPTONIC DECAYS 

The entangled B° — B ° mesons, produced in the decay 
of the T(4S) resonance, can both decay semi-leptonically. 
The effects of decoherence on the resulting dilepton signal 
was studied in Here we calculate these effects using 


the formalism described in the previous section. The 
entangled B d — B ° state can be written as 


\m) 




( 8 ) 


The time evolution of the above state is described by the 
following density matrix 0-0]: 


p(ti,t 2 ) = 0 p2it 2 ) + P2{tl) 0 pi{t 2 ) 

~Pz{ti) 0 Piit 2 ) - Pi{t\) 0 /0 3 (f 2 )) , (9) 


where pi(t) = pBo(t), p 2 (t) = pgo{t) which are given 
in Eq. 0, while p 3 / 4 (t) = Ei -^ 3 / 4 ( 0 )^, where 
P 3 / 4 ( 0 ) = \P°(B 0 )) and are given by 


P3{t) 

kp-T t 
2 e 

Pi(t) 

lp-rt 


-a sh - ie xt a s 
dch a c 

0 

-a sh + ie~ xt a s 

dch + 6 d c 

0 


dch + 6 a c 
—a s h + ie~ xt a s 

0 


dch 6 d c 
—a sh - ie~ xt a s 


0 



° \ 

0 ID) 

2 d s h J 


Here the parameters are as in Eq. 0- The double decay 
rate, G(f, t±;g, t 2 ), that the left-moving meson decays at 
proper time t\ into a final state /, while the right-moving 
meson decays at proper time t 2 into the final state g , 
is then given by Tr [(Of 0 O g ) p(t\, t 2 )\. From this a 
very useful quantity called the single time distribution, 
T(/, ^;t), can be defined as T(f,g;t) = / 0 °° dr G(f,r + 
t',g,r), where t = ti — t 2 is taken to be positive. 

We now consider the decays of B ° mesons into semilep¬ 
tonic states h 1 12 , where h stands for any allowed charged 
hadronic state. Under the assumption of CPT conserva¬ 
tion and no violation of A B = A Q rule, the amplitudes 
for B d /B® into h~l + u can be written as 


A ( B° d -► h~l + v) = M h , A -> h~l + v ) = 0 , 

(U) 

whereas the amplitudes for B®/B d into h + l~v are 

A(B° d ^h + rD)= 0, A -*• /i+rp) = M£. 

( 12 ) 

There are two important observables which can be af¬ 
fected by interaction with the environment. One is the 
ratio of the total same-sign to opposite-sign semileptonic 
rates 


_ r(h+,h+) + r(h-,h-) 
d r(h+,h-) + r(h-,h+) ’ [ } 

and the other is the total B° — B d mixing probability 

r(h+,h + ) + r(h-,h~) 

Xd ~ T(h+, h+) + T{h-,h-) + T(h+,h~) + T(h~,h+) ' 

(14) 
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Time independent probabilities, T(f, g), can be obtained 
by integrating the distribution T(f,g- 1 t) over time. 

The expressions for R d and Xd are obtained to be 


Rd = 


1 - (i- y 2 3 4 ) ((i + A ') 2 +x 2 ) 


-1 


-1 ’ 


1 + (1 — y 2 ) ((1 + A') 2 + x 2 ) 

Xd = \[i^{i-y 2 ){ii + \') 2 +x 2 ) 


(15) 


(16) 


where we x = Am/T, y = AT/2T and A' = A/T. We 
see that Rd and Xd are both functions of (1 — y 2 ) and 
(1 + A') 2 . It is interesting to note that in the limit of 
small \' and y , these combinations have a linear term in 
A' but only a quadratic term in y. Thus we see that along 
with Aro^ and AT^, these observables also depend upon 
the decoherence parameter A. 

For the observable Rd, the last experimental update 
was given about two decades ago fgj. This value was used 
in ref. Q to estimate the value of A to be (—0.072±0.118) 
ps -1 . It is important to reanalyze the BaBar and Belle 
data on the time dependent mixing asymmetry in terms 
of the three parameters (A, Amj, AT^) using the expres¬ 
sion given in Eq. ©• One should also obtain the value 
of Xd from CDF, DO and LHCb. Then the expression in 
Eq. (ITGT) can be verified using the values obtained from 
the fit to the time dependent mixing asymmetry. This 
will provide an additional consistency check on assump¬ 
tions made regarding decoherence. Finally, the values of 
A, A rrid and AT^ from the A m i x (i, A) fit can be used in 
Eq. (J3J to obtain a clean measurement of sin 2/3. 

The present analysis can easily be extended to the 13° 
system as well. The expression for the time dependent 
CP asymmetry in the mode 13° —> J/if^ will be a func¬ 
tion of four parameters: A, sin2/3 s , A m s and Ar s . The 
time dependent mixing asymmetry defined in Eq. ([7]) will 
determine A, Am s and AT S . These two time-dependent 
asymmetries should be re-analysed using a four param¬ 
eter fit for a clean determination of sin2/3 s , Am,, Ar s 
and A. Also, like sin 2/3 d, the extraction of sin 2[3 S from 
time dependent CP asymmetry in the mode 13° —> J/il>4> 
is restricted due to penguin pollution. In this case, the 
analysis of CP violation is more involved in comparison 
to 13° —> J/ifKs- This is due to the fact that the fi¬ 
nal state involves two vector mesons. The admixture of 
different CP eigenstates can be disentangled through a 
time-dependent angular analysis of the decay products of 
the vector mesons |3ll . HH. The penguin contribution to 


B ° —> J/ibd> can be estimated using decays —> J/ibp 

a,nd B° U J/iPK* \mm- 


V. ESTIMATION OF A: AN EXAMPLE 

Here we make an attempt of a clean determination of 
A, A rrid and Ar^ using the experimental data of the time 
dependent flavor asymmetry of semi-leptonic B° decays 
as given in Ref. 0. We perform a x 2 fit to A m j x (At, A), 
using the efficiency corrected distributions given in Ta¬ 
ble I of Ref. [qJ. First, the fit is done by assuming no 
decoherence, i.e., A = 0. In this case, we find Amj = 
(0.489T0.010) ps" 1 and Ar d = (0.087±0.054) ps” 1 with 
X 2 /d.o.f = 8.42/9. We then redo the fit including deco¬ 
herence. This gives A — (—0.012±0.019) ps -1 along with 
A m d = (0.490 ± 0.010) ps" 1 and Ar d = (0.144 ± 0.088) 
ps -1 with x 2 /d.o.f = 8.02/8. Thus we see that the deco¬ 
herence parameter A is very loosely bounded. The upper 
limit on A is 0.03 ps -1 at 95% C.L. We also find in this 
example that A m d is numerically unaffected where as 
Ar^ can be affected by inclusion of decoherence. Given 
the wealth of data coming from LHCb and expected from 
the KEK Super B factory, a clear picture is expected to 
emerge. 


VI. CONCLUSIONS 

In this work, we have studied the effect of decoherence 
on two important observables sin 2/3 and A m d in a neu¬ 
tral meson system. We find that the asymmetries which 
determine these quantities are also functions of the deco¬ 
herence parameter A. Hence it is imperative to measure 
A for a clean determination of these quantities. We sug¬ 
gest a re-analysis of the data on the above asymmetries 
for an accurate measurement of all the three quantities 
A, sin 2/3 and A m d - The present analysis can easily be 
extended to the 33° system as well. 

Acknowledgments .— We thank Kajari Mazumdar for 
helpful discussions on several parts of this analysis. We 
thank David Hitlin for suggesting an example for deter¬ 
mination of A. We also thank Joaquim Matias, David 
London and B. Ananthanarayan for useful comments. 
The work of AKA and SB is supported by CSIR, Gov¬ 
ernment of India, grant no: 03(1255)/12/EMR-II. 


[1] S. W. Hawking, Commun. Math. Phys. 87 (1982) 395. 

[2] J. R. Ellis, J. S. Hagelin, D. V. Nanopoulos and M. Sred- 
nicki, Nucl. Phys. B 241, 381 (1984). 

[3] U. Weiss, Quantum Dissipative Systems, Third Edition 
(World Scientific 2008). 

[4] H.-P. Breuer and F. Petruccione, The Theory of Open 
Quantum Systems (Oxford University Press 2002). 


[5] S. Banerjee and R. Ghosh, Phys. Rev. E 67, 056120 
(2003). 

[6] S. Banerjee, A. K. Alok and R. MacKenzie, 
arXiv: 1409.1034 [hep-ph], 

[7] R. A. Bertlmann and W. Grimus, Phys. Rev. D 64, 
056004 (2001) hep-ph/0101160 . 

[8] H. Albrecht et al. [ARGUS Collaboration], Phys. Lett. B 
324, 249 (1994); J. E. Bartelt et al. [CLEO Collabora- 






5 


tion], Phys. Rev. Lett. 71, 1680 (1993). 

[9] A. Go et al. [Belle Collaboration], Phys. Rev. Lett. 99, 
131802 (2007) quant-ph/0702267 [QUANT-PH]]. 

[10] G. Lindblad, Commun. Math. Phys. 48, 119 (1976). 

[11] V. Gorini, A. Kossakowski and E. C. G. Sudarshan, J. 
Math. Phys. 17, 821 (1976). 

[12] R. Alicki and K. Lendi, Quantum Dynamical Semigroups 
and Applications (Lect. Notes Phys. 717 (Springer, Berlin 
Heidelberg 2007)). 

[13] E. C. G. Sudarshan, P. M. Mathews and J. Rau, Phys. 
Rev. 121, 920 (1961); K. Kraus, States, Effects and 
Operations: Fundamental Notions of Quantum Theory 
(Springer Verlag 1983). 

[14] A. K. Alok et al, Work in progress. 

[15] R. Fleischer, Phys. Rev. D 60, 073008 (1999) 

hep-ph/9903540 . 

[16] K. De Bruyn and R. Fleischer, JHEP 1503, 145 (2015) 
arXiv: 1412.6834 [hep-ph]]. 

[17] R. Fleischer, Eur. Phys. J. C 10, 299 (1999) 

hep-ph/9903455 . 

[18] K. A. Olive et al. [Particle Data Group Collaboration], 
Chin. Phys. C 38, 090001 (2014). 

[19] G. Abbiendi et al. [OPAL Collaboration], Phys. Lett. B 
493, 266 (2000) hep-ex/0010013 . 

[20] D. Buskulic et al. [ALEPH Collaboration], Z. Phys. C 
75, 397 (1997). 

[21] J. Abdallah et al. [DELPHI Collaboration], Eur. Phys. 
J. C 28, 155 (2003) liep-ex/0303032 . 

[22] M. Acciarri et al. [L3 Collaboration], Eur. Phys. J. C 5, 
195 (1998). 


[23] F. Abe et al. [CDF Collaboration], Phys. Rev. D 60, 
072003 (1999) hep-ex/9903011. 

[24] B. Aubert et al. [BaBar Collaboration], Phys. Rev. D 73, 
012004 (2006) hep-ex/0507054 . 

[25] K. Abe et al. [BELLE Collaboration], Phys. Rev. D 
71, 072003 (2005) [Erratum-ibid. D 71, 079903 (2005)] 
hep-ex/0408111 . 

[26] V. M. Abazov et al. [DO Collaboration], Phys. Rev. D 
74, 112002 (2006) liep-ex/0609034 . 

[27] R. Aaij et al. [LHCb Collaboration], Phys. Lett. B 
709, 177 (2012) arXiv:1112.4311 [hep-ex]]; R. Aaij et 
al. [LHCb Collaboration], Eur. Phys. J. C 73, no. 12, 
2655 (2013) arXiv:1308.1302 [hep-ex]]; R. Aaij et al. 
[LHCb Collaboration], Phys. Lett. B 719, 318 (2013) 
arXiv:1210.6750 [hep-ex]]. 

[28] P. Huet and M. E. Peskin, Nucl. Phys. B 434, 3 (1995) 
hep-ph/9403257 . 

[29] F. Benatti and R. Floreanini, Phys. Lett. B 465, 260 
(1999) hep-ph/9909361 . 

[30] F. Benatti, R. Floreanini and R. Romano, Nucl. Phys. B 
602, 541 (2001) hep-ph/0103239 . 

[31] A. S. Dighe, I. Dunietz, H. J. Lipkin and J. L. Rosner, 
Phys. Lett. B 369, 144 (1996) hep-ph/9511363 . 

[32] A. S. Dighe, I. Dunietz and R. Fleischer, Eur. Phys. J. C 
6, 647 (1999) hep-ph/9804253 . 

[33] S. Faller, R. Fleischer and T. Mannel, Phys. Rev. D 79, 
014005 (2009) arXiv:0810.4248 [hep-ph]]. 




